Pseudo-gap and vertex correction of electron-phonon interaction 
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The strong-coupling Eliashberg theory plus vertex correction is used to calculate the maps of 
transition temperature (T c ) in parameter-space characterizing superconductivity. Based on these T c 
maps, complex crossover behaviors are found when electron-phonon interaction increases from weak- 
coupling region to strong-coupling region. The doping-dependent T c of cuprate superconductors and 
most importantly the pseudo-gap can be explained as the effects of vertex correction. 
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The standard strong-coupling theory has no bound 
on T c . Recently, the Eliashberg functions a 2 F(ui) ex- 
tracted from the measurements of infrared optical con- 
ductivity [l[ and ARPES spectrum for copper-oxides 
superconductors predicted very strong electron-phonon 
interaction and very high T c over the experimental val- 
ues [l| . The T c in mean-field approximation of Eliashberg 
theory is higher than experimental T c . In the situation 
of very strong electron-phonon coupling, non-adiabatic 
effects of electron-ion system will be so important that 
the electrons are dressed heavily by lattice vibrations and 
the conventional strong-coupling theory needs to be gen- 
eralized to include the non-adiabatic effects or the ver- 
tex corrects beyond Migdal's theorem [3j-|5[. The be- 
havior of crossovers when electron-phonon interaction 
evolving from weak-coupling region to strong-coupling 
region were found in different theoretical calculations [g- 
8]. These crossovers are expected to prevent T c from 
infinitely increasing with electron-phonon interaction. If 
the electron-phonon interaction is the underlying pairing- 
mechanism for unconventional superconductors such as 
the cuprate superconductors, it should provide reason- 
able explanations of pseudo-gap and the dome-shape of 
doping-dependent T c . 

In previous works [1, 0, [l(|, we predicted that the 
highest T c of cuprate superconductors is close to 160K 
consistent with the present record of cuprate supercon- 
ductors [f| and successfully explained the spatial anti- 
correlation between energy gap and phonon energy for 
cuprate superconductor Bi2212 The up-limit of T c 
for iron-based superconductors is about 90K [j| and there 
are 34-35 K space left to increase T c beyond present 
record about 55-56 K [ll[. If the effects of vertex cor- 
rection (or non-adiabatic effects) are strong, it's very 
difficult to realize home-temperature superconductor in 
high-pressure metal hydrogen and other hydrogen-rich 
materials [10| | . 

In this paper, the T c maps including the influences 
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of vertex corrections are studied. Complex crossovers 
are found on these T c maps when the parameter A of 
electron-phonon coupling increases from weak-coupling 
region to strong-coupling region. These crossovers are 
very close to the well known A=2 at which the value 
of T c reaches its maximum [13]. The pseudo-gap and 
phase diagram with same topology as the phase diagram 
of doping-dependent T c for cuprate superconductor are 
explained as the effects of vertex correction. The inter- 
play of vertex correction and Coulomb interaction can 
suppress theoretical T c to access experimental values [l[ 
The calculations of vertex corrections are greatly 
simplified under isotropic approximation because the 
electron-phonon interactions are included in the vertex 
corrections only by the functions of electron-phonon in- 
teraction A n defined as A„ = 2 J Q dwa 2 F(oj)uj / '{uj 2 + 
cj^). When temperature is very close to T c 
the energy-gap equation 0, [|| is simplified to 
Sn'=-oo K nn '(A n > /\u! n '\) = 0. The kernel matrix is ex- 
pressed as 

■Knn' — [Art— n'^nn' ~f" C nn '^d n ' S nn ' H n ' , (1) 

TT _ V"^ f S n ' n "\uJn"\ . . , 

fin' / I i m ' A n i — n " Si n ' n " S n ' S n " & n " \ . 

z — ' TTKbI 

n" — — oo 

where the parameters A nn i = 1 - V r ^ n , , B n „> = 1 - V^ n , , 
s n = w n /\w n \ and a n = (2/w) arctan(i? B /Z n |w n |). The 
three parameters of vertex correction V^ n ,, V^ n , and 
C nn i can be found in Ref.[f|. The Coulomb pseudo- 
potential is defined as /i* = ^q/(1+iiq\jx[Eb/u)q))-, where 
Ho = N(0)U, U the Coulomb interaction between elec- 
trons and wo characteristic energy of typical phonon cor- 
related to superconductivity. If the vertex corrections 
are ignored, three parameters V^ n ,, V^ n , and C nn ' are all 
equal to zero and the kernel Eq.([T|) of energy-gap equa- 
tion reduces to the general form without vertex correc- 
tion after some symmetrizations and simplifications. 
It's convenient that the K nn i matrix is symmetrized as 
in Ref. [14|. The Eliashberg functions a 2 F(u>) have the 
same approximation as in Ref.[12j. Other details in our 
calculations can be found in Ref[5|, [l(| 

The parameter ftp/ Eb measures and controls the mag- 
nitude of vertex correction in perturbing calculation. 
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FIG. 1: The evolution of T c map on A-f2p plane with in- 
creasing strengths of vertex corrections (decreasing effective 
band-width E B ) with (a) E B = oo, (b) E B = 1.7 eV and (c) 
Es = l eV. The Coulomb pseudo-potential fi* =0.1. 



From electron point of view, the vertex correction or non- 
adiabatic effect can be controlled by the effective band- 
width Ep, on the other hand, from ion point of view, it 
can be controlled by the cutoff ujq of phonon energy or 
£!p in Einstein model. In this work, the vertex correc- 
tion is controlled by the effective band-width Ep within 
the range from 0.5 eV to 5 eV. The situation Ep = oo 
is equivalent to no vertex correction. The smaller Ep 
means possible the stronger vertex correction. 

The FigQJa,b,c) illustrate the evolution of T c map on 
X-ilp plane with decreasing Ep. The FigOJa) is the T c 
map having been obtained in the previous work without 
considering vertex corrections Q. When Ep=1.7 eV, the 
large deformation of T c map with strong vertex correc- 
tion is shown in FigUJb) near the well known A=2.0 in 
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FIG. 2: (a) The T c change along two arrows shown in Fig[T|c) 
with fixed phonon energies f2p=80 meV and 95 meV respec- 
tively, (b) The T c change along two arrows shown in Fig[3jc) 
with fixed Coulomb pseudo-potentials fi*=0.10 and 0.25 re- 
spectively. 
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FIG. 3: The evolutions of T c map on fj*-\ plane (f2p=69 
meV) with decreasing effective band- width (a) E B = oo, (b) 
E s = 1.7 eV and (c) E s = 1.0 eV. 



the region of high phonon energy. With Ep decreasing 
to 1 eV further, the region with strong vertex correc- 
tion rapidly expands and occupies large part of param- 
eter space with fip>80 meV in FigJIJc). In the region 
fip<80 meV, the T c is strongly suppressed however there 
are no discontinuous changes of T c or breaking of con- 
tour lines. An important result from the Fig[T]is that T c 
does not change with A monotonously if phonon energy 
fip is high enough. The FigJ5Ja) shows the changes of 
T c with A along two arrows A and B shown in FigJIJc). 
If Slp=80 meV, the T c monotonously increases with A. 
However for fip =90 meV, the T c first increases with 
A, reaches the maximum at A ~ 1.5 — 1.7 and then 
quickly decreases with increasing A. Further increasing 
A>2, T c will be very low due to strong vertex correc- 
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tions. The non-monotonous A-dependent T c in Figj2ja) 
had been found in the non-adiabatic theory of supercon- 
ductivity !(j . Some crossover behaviors from weak cou- 
pling to strong coupling region had been predicted in 
Holstein-Hubbard model solved numerically by quantum 
Monte Carlo method 0] and in polaron theory Q. It's 
very reasonable that the non-monotonous A-dependent 
T c is equivalent to the crossovers found in QMC calcu- 
lation [3] and polaron theory Q. So only the leading 
vertex correction can describe qualitatively very well the 
electron-phonon interaction in strong coupling region. 
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FIG. 4: (a) The T c map on E B -X plane with n*=0.25 and 
Q,p=72 meV. (b) The open circle line is the evolution of T c 
from PI to P2 in (a) but fi* linearly decreases from 0.3 to 
0.1. The solid line T* is the standard results in strong cou- 
pling theory without vertex correction, (c) The S-X relation 
is adopted in Ref.[lJ 

The Fig|3][a) is the normal T c map on /P-A plane with- 
out vertex correction |5j. The figure shows that when 
/P>0.2, T c is insensitive to the change of /P. The break- 
ing contour lines with T c =0 K are because of the inac- 
curate calculations when T c <0.1 K if only iV=200 Mat- 
subara energies are used. The contour lines with T c >0.1 
K are accurate enough. If the Coulomb pseudo-potential 
and vertex correction work together, the situation will 



change drastically and some new interesting results will 
appear. The large deformations arc found in FigJ3Jc) if 
Eb decreases to 1.0 eV. As expected, the large deforma- 
tions and discontinuous changes of contour lines appear 
on the T c map when /P>0.20. The contour lines with iso- 
values from T c =20 K to 200 K are packed together within 
the rectangle region in FigEtc) with 0.15</P<0.25 and 
A>2. The figure clearly shows that if the Coulomb 
pseudo-potential /P is larger enough, the T c will change 
with A non-monotonously. The changes of T c along two 
arrows with ^t*=0.1 and 0.25 are plotted in FigJ^b). For 
/P=0.25, T c first increases with A until reaches the max- 
imum at A=2.2 and then sharply decreases to smaller 
value at A=2.5. The crossover behavior is enhanced by 
strong Coulomb interaction. 
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FIG. 5: (a) The effects of three non-adiabatic parameters 
V^ n i y„ n ' and C nn i on T c . (b,c) The changes of average values 
of diagonal and off-diagonal matrix elements of \A nn i \ = |1 — 
v nn'l \ B nn' \ = I 1 ~ V^n' I and \ G nn> I with increasing A. 

The T c map on E^-A plane is presented in FigHta) 
with fip=72 meV. If E^ increases but A keeps un- 
changed, the T c monotonously increases with E# until 
to the limit of non- vertex-correction. More interestingly, 
on this map, the T c is non-monotonous dependent on 
E# along straight line from PI to P2 companying by 
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the decrease of A from 5.0 to 0.2. The non-monotonous 
dependence of T c on effective band-width is equiv- 
alent to the band- filling effect of T c . Our results show 
that, if fip >80 meV, the suppression of Tc will be more 
prominently. 

The values of T c obtained from standard strong- 
coupling theory are generally higher than those mea- 
sured in experiments. The copper-oxides superconduc- 
tors BioS^CaC^Os+a and Bi2Sr2Cu20g+5 studied in 
Ref [H [If have very strong electron-phonon interactions 
A~2.36-2.85 and overestimated T c in underdoped sam- 
ples. With increasing doping S, the values of A decrease 
to about 0.35-1.42 [lj. The effective band-widths E B of 
conducting electrons for these cuprates are distributed 
from 1 eV to 3 eV. The effective phonon energies are dis- 
tributed from 50 meV to 80 meV. We re-calculate the 
values of T c along straight line P1-P2 in FigfJJa) under 
assumption that the Coulomb interaction is strong in un- 
derdoped region /i*=0.3 at PI and weak in overdoped re- 
gion /i*=0.1 at P2. For simplicity, /x* linearly decreases 
from 0.3 at PI to 0.1 at P2. As shown in Figgjb), if 
A <4.0, the values of T c are reduced from around 200 K 
to lower than 150 K and close to experimental values [lj]. 
In strong-coupling region 4.0> A >3.0, T c is very low. 
Our results are provided an explanation to pseudo-gap in 
underdoped region shown in Figf^b). The cooper-pairs 
pre-form at T* the transition temperature of mean field 
theory (MFT) or the standard strong-coupling theory. 
However strong non-adiabatic effects induce the instabil- 
ity of cooper-pairs and the real T c has lower value. The 
T* degenerating with T c in overdoped region is similar 
to the example (1) of Fig.9 in Ref.[15[. The pre-formed 
cooper-pairs in cuprate superconductors are supported 
by measurements of Nernst effect [l6[ , specific- heat [TtJ 
and many other methods. 

An interesting result is that at very strong coupling 
A >4.0, the effects of vertex corrections superficially 
become weak. Even there are positive vertex correc- 
tion that had been found in other work The elec- 
tronic states in region SP with A >4.0 are strong- 



coupling pairs [Tg ■ The FigfJJb) shows a crossover from 
BCS state to strong-coupling pairs state with increas- 
ing electron-phonon interaction A. It's obviously that 
FigHfb) has the same topology as the well-known T c - 
doping phase diagram. It's dependent on whether the 
parameter A electron-phonon interaction decreases with 
increasing doping or not. This point had been proved in 
recent experiments [l|, 0. The S — A curve in FigfJIc) 
is based on data in Ref.[l[. It's urged that there will be 
other experiments supporting this point. 

In order to analyze our results more deeply, we present 
individually the effects of non-adiabatic parameters V^ n , , 
V£ n , and C nn > on T c in FigUfa). The T c -A curve labeled 
with V B is calculated by allowing V^ n , ^ and setting 
V^ n ,=0 and C nn .=Q. Other curves are obtained with the 
same manner. We can find that the dome shape curve 
of T c in the region A <4.0 is generated by the effects of 
V^ n , . In the region A >4.0, the effects of V^ n , and C nn > 
cancel the effects of V^ n , so that the strong coupling pairs 
in mean field approximation are restored and even have 
higher T c . This fact can be clarified from the FigJSJb,c) 
in that the averages of the absolute values of diagonal and 
off-diagonal elements of parameter matrix |A„„/ = |1 — 
Vnn'l \ B nw\ = |1 - Km' | and \C nn > \ have larger changes 
mainly in the region 2.5< A <4.0. Moreover in the region 
A >4.3, the I I and |-B n ra'| are close to normal values 
1.0 just as in weak coupling region. Additionally, the 
average values of diagonal elements of parameter matrix 
\C nn i I and \B nn i\ steadily increase with A and lead to 
positive vertex-correction. 

In summary the non-monotonous changes of T c with 
increasing A show the crossover behaviors near A=2 when 
A evolving from weak-coupling region to strong-coupling 
region. The crossovers can explain both the pseudo-gap 
phenomenon and the dome shape of doping dependent T c 
of cuprate superconductors. The T c maps in the previ- 
ous paper [f| and the maps with vertex corrections in this 
paper provide very comprehensive understanding of su- 
perconductivity of superconductors. The author thanks 
Prof. E. Cappelluti for very helpful discussions. 
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